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ABSTRACT 

Using the Australia Telescope Compact Array, we have carried out a survey of the H I emission in the direc- 
tion of the "barrel-shaped" supernova remnant (SNR) G320.4-1.2 (MSH 15-52) and its associated young pulsar 
B1509-58. The angular resolution of the data is 4'0 x 2.'7, and the rms noise of the order of 30 mJy/beam (~ 0.5 
K). The Hi observations indicate that the N-NW radio limb has encountered a dense Hi filament (density ~ 12 
cm -3 ) at the same LSR velocity than that of the SNR (Vlsr ~ -68 km s" 1 ). This H I concentration would be 
responsible for the flattened shape of the NW lobe of G320.4-1.2, and for the formation of the radio/optical/X- 
ray nebula RCW 89. The emission associated with the bright knots in the interior of RCW 89 can be explained 
as arising from the interaction between the collimated relativistic outflow from the pulsar and the denser part of 
this Hi filament (density ~ 15 cm" 3 ). The S-SE half of the SNR, on the other hand, seems to have rapidly 
expanded across a lower density enviroment (density ~ 0.4 cm" 3 ). The H I data also reveal an unusual H I feature 
aligned with a collimated outflow generated by the pulsar, suggestive of association with the SNR. The anomalous 
kinematical velocity of this feature (Vlsr ~ 15 km s" 1 ), however, is difficult to explain. 

Subject headings: ISM: individual (G320.4-1.2, RCW 89, MSH15-52) — ISM: structure — pulsars: individual 
(B1509-58) — radio lines: ISM — supernova remnants 
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1. INTRODUCTION 

The radio morphologies of supernova remnants (SNRs) are 
usually approximated by limb-brightened spherical shells, in 
which the shock front accelerates electrons to synchrotron- 
emitting energies. However, observations at increasingly higher 
resolution and sensitivity have demonstrated that actual SNR 
morphologies are much more complex, exhibitin g "break-out", 
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■ ico n 1996|; |Gaensler 199"8f Dubner et al. 
1998b). Such a diversity reflects not only variations in the prop- 
erties of the progenitor star and of the explosion itself, but also 
echoes the physical conditions of circumstellar and interstellar 
gas and magnetic fields. The interpretation of SNR morpholo- 
gies has the additional complication that the images observed 
are a two-dimensional projection of a three-dimensional object. 
Radio spectroscopy of the interstellar medium (ISM) provides 
a valuable tool to unravel these 3D effects using kinematic sig- 
natures. As part of an on-going e ffort to study the inter action 
of SNRs with t heir environment ([Dubner et al. 19984 D ub- 
ner et al. 1998b; Dubner et al. 1999| ; |Giacani et al 2000| ), we 
present a study of the distribution and kinematics of the neutral 
hydrogen in the vicinity of SNR G320.4-1.2 (MSH 15-52). 

SNR G320.4-1.2 has an unusual radio morphology, consist- 
ing of two loosely connected emission regions, one to the N- 
NW and the other to the S-SE. Figure |l] shows the radio and 
X-ray emission from the SNR. The young pulsar PSR B1509- 
58 lies between the two main radio components of the SNR; 
the pulsar's position as measured by Gaensler et al. (1999), 
hereafter G99, is marked in Figure [j] with a cross. The 



optical/radio/X-ray nebula RCW 89 is located on the NW ex- 
treme of G320.4-1. 2. 

Like many othe r SNRs, G320.4-1.2 has a " b ilateral" or "bar - 
rel" morphology ( Kesteven & Caswell 1987| ; Uaensler 1998). 
Many different mechanisms have been proposed to account for 
this morphology, relating to both th e progenitor star and super- 



nova explosion (" i ntrinsic effects") ([Kesteven & Caswell 1987 
Manchester 1987|;[Brighenti & D'Ercole 1994|; W illingale"etal 



1996) and to subsequent interaction with the SNR environment 
("ext rinsic effects") (Bisnovatyi-Ko gan, Lozinskaya, & Silich 
1990; |Fulbright & Reynolds 1990[ ). Most recently, Gaensler 
(1998) has argued that the bilateral appearance arises when a 
SNR expands into an elongated low-density cavity; the radio 
emission then originates at the sites where the shock encoun- 
ters the walls of the cavity. 

This peculiar SNR and its energetic young pulsar have been 
extensively observed throughout the electromagnetic spectrum 
(see G99 for a summary of previous observations). In par- 
ticular, recent X-ray and radio data provide convincing evi- 
dence that the SNR and pulsar resulted from a single super- 
nova explosion which occurred approximatel y 1700 yrs ago 
( [Tamura et al 1996|; |Brazier & Becker 1997|; G99; G aensler 
et al. 2001). Specifically, argue that the entire complex 
can be interpreted as the result of a low-mass or high-energy 
explosion occurring near one edge of an elongated low density 
(n ~ 0.01 cm" 3 ) cavity. To the NW the SNR has expanded into 
dense material (n ~ 1 -5 cm" 3 ), while to the SE the remnant has 
expanded rapidly across the low density cavity. Meanwhile, the 
pulsar appears to be generating twin collimated outflows, the 
northern of which is interacting with the SNR to produce a col- 
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lection of compact radio/X-ray knots embedded in RCW 89. A 
distance of (5.2 ± 1 .4) kpc was derived for G320. 1-1 .2 by G99, 
based on H I absorption measurements, in agreement with the 
results of Caswell et al. (1975). 

Within this complex scenario, it is clear that the interpreta- 
tion of the existing body of observational data can greatly bene- 
fit from an accurate knowledge of the distribution and kinemat- 
ics of the surrounding neutral hydrogen emission. Here we re- 
port on an H I survey of SNR G320.4-1.2 and its environment, 
carried out with the Australia Telescope Compact Array. 

2. OBSERVATIONS AND DATA REDUCTION 

Interferometric Hi observations towards G320.4-1.2 were 
carried out using the Australia Telesco pe Compact Array 
(ATCA; prater, Brooks, & Whiteoak 1992 ), a six-element syn- 
thesis telescope located near Narrabri, NSW, Australia, during 
a session of 13 hours on 1998 October 13. In order to achieve 
optimal u—v coverage in the radial direction, the array was used 
in the non-standard 210-m configuration, in which five anten- 
nas are positioned closely together so as to provide 10 baselines 
ranging from 31 to 214 m. The sixth (fixed) antenna is sepa- 
rated from the other five elements by 5-6 km, and was not used 
in our analysis. 

A region surrounding the SNR was surveyed in a mosaic of 
19 pointings, following a hexagonal grid to cover six square 
degrees . The separation between grid-points of about 16' satis- 
fies the Nyquist sampling criterion for the primary beam width 
of 33' (FWHM). Observations were centered at a frequency 
of 1420 MHz, using 1024 channels over a total bandwidth of 
4 MHz. The absolute flux density scale was determined by 
assuming a flux density for PKS B 1934-638 of 14.85 Jy at 
1.42 GHz. Observations of PKS B1934-638 were also used 
to determine the spectral response across the observing band. 
Time variability in the antenna gains was calibrated using regu- 
lar observations of PKS B 1540-828. After editing and calibra- 
tion, continuum emission was subtracted from the data in the 
u — v plane u sing the task UVLSF in AIPS (v an Langevelde & 
Cotton 1 990). The data were th en transferred to the MI RI AD 
package (Sault & Killeen 1999), where all pointings were im- 
aged simultaneously to form 390 mosaicked planes, covering 
LSR velocities between -150 and +171 km s" 1 at a velocity 
resolution of 1 km s" 1 (channel separation of 0.82 km s" 1 ). 
This cube was then deconvolved using a mosaicked maxim um- 
entropy technique (Sault, Staveley-Smith, & Brouw 1996), and 



then convolved with a gaussian restoring beam. 

Any analysis of H I emission requires measurements at all 
spatial scales down to the resolution limit. However, our ATCA 
observations are not sensitive to structures larger than 15'-20'. 
To recover this short-spacing data, we have utilized H I obser- 
vations made using the Parkes 64-m radio-telescope, carried out 
as part of the Southe rn Galactic Plane Survey (SGPS; M cClure- 
Griffiths et al. 2001) with an angular resolution of 14'. The 
Parkes data were multiplied by a factor of 0.9 to match the 
flux calibration applied to the ATCA data; interferometric and 
single-dish data were then appropriately weighted according 
to their respective primary beam shapes and combin ed using 
the task IM MERGE in MIRIAD ([Stanimirovic 1999|; Sault & 
Killeen 1999). The final combined cube has units of Jy beam" 1 . 
To convert to units of brightness temperature, these data must 
be multiplied by a factor of 16.3 K Jy" 1 beam. The angular res- 
olution of the data is 4.'0 x 2.'7 (Position Angle 55°), and the 
rms in line-free channels is ~ 30 mJy beam" 1 (~ 0.5 K). 



3. RESULTS 
3.1. General H I distribution 

Figure ^ illustrates the general properties of H I emission in 
the direction of G320.4-1.2. The top panel shows an aver- 
age H I profile of the observed region taken from the Parkes 
SGPS data, while the bottom panel displays the circular rota- 
tion model towards I = 320°, b = — 1°, assuming Rq = 8.5 kpc 
and 9o = 220 km s" 1 and using the rotation curve of Fich, Blitz 
& Stark (1989). 

The neutral hydrogen emission in this direction has a com- 
plex structure, and is spread over more than 200 km s" 1 . The 
five main peaks of this spectrum can be interpreted as emission 
from gas related to the Sagitarius-Carina galactic arm (that the 
line of sight crosses at both Vlsr ~ -50 and ~ +80 km s" 1 ), the 
Scutum -Crux arm (at velocities of Vlsr ~ _ 20 and + 40 km 
s" 1 ), and to the local gas contribution near km s" 1 . The SNR 
would be located on the far border of Carina arm (Georgelin et 
al. 1987). 

We investigate the H I emission distribution in the entire ob- 
served velocity interval, looking for morphological and/or kine- 
matical evidences of interaction of G320.4-1.2 and its different 
components with the surrounding H I gas. Figure |] provides 
an overview of the distribution of atomic hydrogen in the sur- 
veyed region. It displays the H I emission between ~ — 8 1 km 
s" 1 and ~ +92 km s" 1 , the interval in which significant H I 
emission is found. Each image is the average of six consecutive 
spectral channels, spanning 5 km s" 1 . However the analysis is 
carried out with the highest possible velocity resolution. The 
field shown is slightly reduced in size with respect to the ob- 
served field to avoid edge effects. The velocity shown in the 
lower left corner of each image corresponds to the velocity of 
the first integrated channel. Superimposed on the H I images 
are a few representative contours of the radio continuum emis- 
sion of G320.4-1.2 at 1.4 GHz as taken from G99, at a spatial 
resolution of 24" x 21". 

From ~ -76 km s" 1 to ~ -66 km s" 1 , the most conspicuous 
feature in H I is an elongated filament seen partially overlapping 
the bright NW part of the SNR. This filament runs diagonally 
from NE to SW, in a direction parallel to the Galactic plane. 
This filament, is discussed in detail in the next section. 

From ~ -60 to ~ -40 km s" 1 , a strong gradient in the gas 
density (assumed to be optically thin) between the NW half of 
the field (closer to the Galactic plane) and the SE corner (away 
from the Galactic plane) is evident. This density gradient per- 
sists up to V ~ -3 1 km s" 1 , although with less contrast than 
at more negative velocities. Significant absorption is evident 
against bright continuum emission from RCW 89 in this veloc- 
ity range (See also G99). 

From ~ -26 to ~ km s" 1 , the H I emission appears quite 
uniformly distributed, with a few brighter concentrations near 
the NW corner. Strong absorption against RCW 89 is observed 
at these velocities. In the most severe cases, this absorption 
is surrounded by residual sidelobes which were not completely 
removed in the deconvolution process. 

Positive velocities along this line-of-sight correspond to dis- 
tant gas, located beyond ~ 12 kpc (see Fig ||). We have exam- 
ined the H I emission at positive velocities, bearing in mind that 
these data can also include contributions from closer gas which 
might be associated with G320.4-1.2 but which has been kine- 
matically disturbed. 

In the velocity range between -1.5 and 8 km s" 1 (i.e. in the 
images labeled -1.5 km s" 1 and +3.5 km s" 1 in Fig. 3), two par- 
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allel bands of H I emission cross over the field from NE to SW. 
The H I concentration that runs close to the SE lobe of G320.4- 
1.2, shows good morphological correspondence with the radio 
continuum emission of the SNR. This H I filament, as seen in 
projection, surrounds the SE lobe of the SNR along the east- 
ern and southern borders. Such correspondence is even closer 
when comparing the H I distribution with the X-ray emission. 
However, since at this kinematical velocity there is no way to 
disentangle the contribution of gas likely to be associated with 
G320.4-1 .2 from the intense local gas emission, we do not con- 
sider this structure as a firm candidate for H I-SNR association. 

From ~ +13 to ~ +22 km s" 1 , there is an unusual, bright 
filament oriented almost perpendicular to the Galactic plane, 
breaking the plane-parallel symmetry observed otherwise. This 
filament overlaps both the brightest region of the SNR and the 
pulsar itself, and will be discussed further in the next section. 

At +33 km s" 1 and +38 km s" 1 , a long, moderately bright 
H I filament is observed in apparent contact with the NW half 
of the SNR. However, since most of the H I in this direction is 
expected to have this plane-parallel distribution, we conclude 
that this feature is most probably distant background material, 
located 14 kpc away. The remaining high- velocity images show 
no features which can be feasibly associated with G320.4-1.2. 

In order to determine the systemic velocity of G320.4-1.2, 
we can use Fig. 3 to analize the characteristics of the HI near 
the SNR. In agreement with previous HI absorption studies ( 
Caswell et al. 1975, G99), the images in Fig. 3 show that HI 
absorption is present towards RCW 89 in the velocity range 
between -61 km s" 1 and -51 km s" 1 . Thus, the kinematical ve- 
locity corresponding to G320.4-1.2 must be more negative than 
-61 km s" 1 . In addition, Hi is observed in emission near the 
northern portion of the remnant up to Vlsr — -76 km s" 1 . We 
therefore conclude that the systemic velocity of G320.4-1.2 lies 
in the range -76 to -61 km s" 1 . Adopting a systemic velocity 
of Vlsr ~ _ 68 km s" 1 , the kinematical distance to the source 
turns out to be about 5 kpc, in good agreement with a previous 
estimate by G99. 

3.2. Particular Hi Features 

In this section we focus on specific H I features, searching 
for favorable H I-SNR associations. 

1. The northern filament in the interval —76 <^ Vlsr <S — 66 km 
s~ l : This velocity interval corresponds to the systemic velocity 
of G320.4-1.2. Figure ^ shows in greyscale the H I distribution 
obtained by integrating between velocities of -76 and -66 km 
s" 1 . The 1.4 GHz radio continuum emission from G320.4-1.2 
(as taken from G99) is plotted in contours. This extended H I 
filament to the NW may represent the dense external wall pre- 
viously suggested by G99 to be slowing down the expansion of 
G320.4-1.2 in this direction. This feature may also be responsi- 
ble for the flattened shape of the NW lobe of the bilateral SNR. 
In particular this filament has a peak in HI emission ( bright- 
ness temperature 80 K at v=-70 km s" 1 ) exactly overlapping 
RCW89. We propose that this H I cloud (about 7' in diameter) 
associated with RCW89, is the dense material with which the 
SNR has colli ded and with which the outflow from the pulsar 
is interacting (Iamura et al. 1996; G99). According to this in- 
terpretation, the bright nebula RCW89 (or, more precisely, the 
bright radio/X-ray knots in its interior) has formed at the site 
of interaction between the collimated relativistic outflow from 
the pulsar and the dense H I cloud. The optical emission associ- 
ated with RCW89 can be explained as arising from the ionized 



outer layers of this high density material, as first suggested by 
Manchester (1987). 

The total H I mass of the extended filament has been esti- 
mated in 4000 M Q , where an appropriate background cor- 
rection has been applied. An average atomic density n~ 12 
cm" 3 has been calculated for this elongated H I feature by as- 
suming that the dimension along the line of sight equals the size 
of the minor axis of the filament. In particular, for the denser 
concentration of this feature, associated with RCW 89, a total 
mass M~ 1700 M Q and a volume density n^ 15 cm" 3 , have 
been derived (in this case a spherical geometry was assumed 
for the H I feature). 

As another evidence of the interaction of G320.4-1.2 with 
its surroundings, we have searched for accelerated shocked H I 
gas all along the northern H I filament, and especially in the di- 
rection of RCW 89, where the major interaction SNR-ISM oc- 
curs. We examined several Hi spectra, as well as RA vs Vlsr 
plots (traced at constant declination), looking for high-velocity 
features, as observed for example in IC 443 (DeNoyer 1977). 
Within the moderate angular resolution of the present data, no 
evidences of accelerated H I gas were found. 

We have also investigated the characteristics of the interstel- 
lar gas where the SE lobe of G320.4-1.2 has expanded. To carry 
out this calculation we have considered a sample volume in a 
location relatively free of H I concentrations around Vlsr ~ ~68 
km s" 1 , the systemic velocity of G320.4-1 .2. The obtained vol- 
ume density is n~ 0.4 cm" 3 . 

The observed contrast in density between the H I near the 
NW lobe and the ambient Hi in the interior of the SNR, is 
consistent with the model of G99 in which the NW half of the 
SNR has encountered the edge of a cavity, while the SE half 
has expanded relatively unimpedded into a lower density envi- 
ronment. 

2. The bright filament perpendicular to the Galactic plane 
in the interval +10 <J Vlsr <, +20 km s" 1 : Figure |] displays in 
detail this H I filament (thick black contours), overlaid on radio 
continuum (greyscale) and X-ray (light grey contours) emission 
from the SNR. The cross indicates the position of PSR B1509- 
58. 

From a morphological point of view, there is a striking align- 
ment of this H I filament with both the NW/SE axis of the X-ray 
nebula and the weak radio features extending south of RCW 89. 
ROSAT and Chandra data demonstrate that the elongated X- 
ray nebula represents opposed collim ated outflows directed (in 



projection) along this NW/SE axis (Brazier & Becker 1997 



Gaensler et al. 2001). G99 argue that the radio emission in 
this region corresponds to a cylindrical sheath surrounding the 
pulsar outflow; it has been shown that the magnetic field in the 
region is well-ordered and ori ented in a direction parallel to the 



axis of the proposed outflow (Milne, Caswell, & Haynes 1993 



G99). It is thus conceivable that a pre-existing column of gas 
with an orientation similar to that of the outflow from the pulsar 
may have provided the conditions necessary to form a cylinder 
and thus generate the observed characteristics. For instance, the 
magnetic field could be frozen to the column of cold H I thus 
explaining its order and orientation. This H I feature may sug- 
gest an association with G320.4-1.2, based solely on the mor- 
phology. However, the high difference between the systemic 
velocity of the SNR and the H I velocity, remains unexplained. 
For this feature a total H I mass M~ 2700 M Q and an average 
atomic density n^ 10 cm" 3 have been estimated. 

4. CONCLUSIONS 
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We have conducted a study of the distribution and kinemat- 
ics of the interstellar neutral hydrogen in a field around the 
bilateral SNR G320.4-1.2. This SNR comprises the bright 
optical/X-ray/radio nebula RCW 89 and the energetic young 
pulsar B1509-58. 

We find several features in H I which suggest that the radio 
morphology of G320.4-1.2 is a result of the distribution of the 
surrounding ISM, combined with the action of the pulsar's rel- 
ativistic outflows. 

The N-NW bright radio limb of G320.4-1.2 shows the ef- 
fects of interaction with an extended H I filament with density 
of ~ 12 cm -3 present at the same velocity as that of the SNR 
(Vlsr ~ _ 68 km s" 1 ). In addition, the interaction of the rela- 
tivistic outflow from PKS B1509-58 with the denser part of this 
filament (n^ 15 cm" 3 ) may have been responsible for the for- 
mation of the bright nebula RCW89 located on the NW extreme 
of G320.4-1.2. The opposite radio limb of the SNR apparently 
expanded into a lower density environment (n ~ 0.4 cm" 3 ) with 
little distortion from a semicircular shape. 

The nature of "barref'-shaped SNRs can be explained as- 
suming that the SNR expands into a low density cavity with 
walls almost parallel, and the bright radio limbs originate where 
the expanding shock encounters the walls (Gaensler 1998). In 
the case of G320.4-1.2, the present H I study has confirmed the 
existence of the H I wall associated with the N-NW lobe. How- 
ever, the observations are not conclusive with respect to the S- 
SE H I counterpart. We have identified H I structures with the 
location and morphology suggestive of an H I-SNR association 
in this direction. However, the fact that they are detected near 
Vlsr ~ +3.5 km s" 1 , indistinguishable in practice from local 
gas contribution, makes the association uncertain. 



Based on morphological coincidences, we suggest a physi- 
cal association between the SNR and an elongated H I feature 
with "Vlsr ~ +15 km s" 1 , oriented perpendicular to the Galac- 
tic plane. This H I structure is strikingly aligned with the pro- 
jected axis inferred for the collimated outflows generated by 
PSR B 1509-58. It is difficult to explain the anomalous kine- 
matical velocity of this feature; we suggest that this gas is in- 
deed at the same distance as is G320.4-1.2, but suffers from 
kinematical distortions. Such perturbations in the velocity field 
are not unexpected in this region, rich in massive stars with 
powerful stellar winds that have stirred the interstellar material 
for thousands of years (Georgelin et al. 1987; Lortet, Georgelin 
& Georgelin 1987). 
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FIG. 1 . — Left: Radio continuum image of the SNR G320.4-1.2 at 1.4 GHz, at a resolution of 24" X 21" (G99). The greyscale varies linearly from to +75 
mJy beam -1 , while contours are at levels of 5, 10, 20, 30, 60, 90 17.0 and 1 50 m ly b eam -1 A few Galactic coordinate lines are included for reference. Right: 
Greyscale and contour ROSAT PSPC X-ray image of G320.4-1.2 ( Trussoni et al. 199( ). The contour levels (in arbitrary units) are 0.5, 1, 1.5, 2, 4 6, 10, 15 and 20. 
The white cross indicates the position of PSR B 1509-58 as determined from mterterometric measurements (G99). 
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FIG. 2. — Top: Average profile of the H I emission as taken from Parkes measurements towards SNR G320.4-1.2, obtained as part of the Southern Galactic Plane 
Survey. Bottom: Galactic rotation curve for / = 320°, b = -1°, according to the Galactic circular rotation model of Fich, Blitz & Stark (1989). 



FIG. 3. — Images of the Hi emission between LSR velocities of ~ —81 and ~ +92 km s , each integrated over six consecutive channels (approximately 5 km 
s _1 ). The velocity shown in the lower left corner of each image corresponds to the first integrated channel; the grey scale for the images varies between +0.08 and 
+6 Jy beam -1 km s . The resolution of the data is 4'.0 X 2'.7 at a position angle of 55° (measured north through east). Superimposed on the H I images are a few 
representative contours of the radio continuum emission of G320.4-1.2 at 1.4 GHz (G99). 
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FIG. 4. — Comparison between radio continuum and Hi in the velocity range -76 to -66 km s . The grey scale, ranging from 15 to 55 mJy beam - ' km 
s -1 corresponds to the Hi distribution , while the contours represent the radio continuum emission at 1.4 GHz (G99.) Plotted contours are 5, 10, 20, 30, 60, 90 
mJy beam -1 . 
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